Superfluidity of "^He Confined in Nano-Porous Media 
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We have examined superfluid properties of *He confined to a nano-porous Gelsil glass that has 
nanopores 2.5 nm in diameter. The pressure-temperature phase diagram was determined by tor- 
sional oscillator, heat capacity and pressure studies. The superfluid transition temperature Tc 
approaches zero at 3.4 MPa, indicating a novel "quantum" superfluid transition. By heat capacity 
measurements, the nonsuperfluid phase adjacent to the superfluid and solid phases is identified to 
be a nanometer-scale, localized Bose condensation state, in which global phase coherence is de- 
stroyed. At high pressures, the superfluid density has a T-linear term, and Tc is proportional to the 
zero-temperature superfluid density. These results strongly suggest that phase fluctuations in the 
superfluid order parameter play a dominant role on the phase diagram and superfluid properties. 

PACS numbers: 67.40.-w 



I. INTRODUCTION 

■^He confined or adsorbed in nanoporous media is an 
interesting model system of strongly correlated Bosons 
under external potential. Recently, we have investigated 
the superfluid and thermodynamic properties of ^He in 
nanoporous Gelsil glass, and have found that the strong 
confinement into the nanopores causes a dramatic change 
in the phase diagraniii^i^iiii^iii^. We show the obtained 
P-T phase diagram in Fig. [1] With increasing pressure, 
the superfluid transition temperature Tc approaches zero 
at 3.4 MPa. Measurements of isochoric pressure have 
suggested that the freezing pressure is at or above 3.4 
MFd^. These behaviors indicate a quantum phase tran- 
sition (QPT) among superfluid, nonsuperfluid and solid 
phases induced by pressure as an external parameter. 

The QPT behavior and the existence of a nonsuper- 
fluid phase between the superfluid and solid phases are 
in striking contrast to the case of bulk "'He. To investigate 
the nature of the nonsuperfluid phase and to understand 
the mechanism of QPT, we have made measurements of 
heat capacity and isochoric pressure'^'^. In this paper, we 
summarize the recent experimental results and propose 
an interpretation that will provide a novel perspective to 
the physics of ^He in porous media: The confinement of 
''He to nanopores fluctuates the phase of superfluid or- 
der parameter, and the phase fluctuations results in the 
localization of Bose-Einstein condensates and a quantum 
phase transition. 
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FIG. f: P — T phase diagram. The superfluid transition tem- 
peratures Tc (green dots) are obtained by torsional oscillatop^ 
studies, and the localized BEG temperatures Tb (red dots) 
and the melting points (red crosses) are by heat capacity^i^. 
Pressure and thermal response measurements^'S, give the melt- 
ing and freezing lines: Tmo; melting onset, Tmt; melting ter- 
mination, Tfo; freezing onset, Tft; freezing termination, re- 
spectively. Arrow indicates the critical pressure Pc at which 
Tc tends to K. 



II. RESULTS: PHASE DIAGRAM AND 
SUPERFLUID PROPERTIES 

A. Phase Diagram 

Here we summarize the results of the measurements 
of pressure and heat capacity, and describe in more 
detail about the torsional oscillator studies, focusing 
on the temperature dependence of the superfluid frac- 
tion. The details of the results have been described 



elsewherei^^^^i^. 

We have employed a porous Gelsil glass^, which is 
manufactured by the sol-gel process. Gelsil has three- 
dimensionally (3D) interconnected nanopores, similarly 
to Vycor. The nominal pore diameter of our glass sam- 
ples is 2.5 nm. Since various pore sizes are available, 
Gelsil has been recently used in helium studies. The 
controllability and wide variety of the pore size were not 
available in Vycor. It was first employed for *He study 
by Miyamoto and Takano (MT)!". They found that the 
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superfluid transition in a 2.5-nm Gelsil sample was de- 
pressed to 0.9 K at ambient pressure. 

We have constructed a heat capacity cell having a ca- 
pacitance pressure gauge^i^. This cell enables us to mea- 
sure the pressure and heat capacity for the same glass 
sample. The sample cell contains four Gelsil disk sam- 
ples (5.5 mm diameter, 2.3 mm thick) which are taken 
from the same batch as the one used in the torsional 
oscillator experiments. 

In the pressure study2,i^^ we measure pressure P{T) 
along isochores. The rates of cooling and warming of the 
cell are also recorded simultaneously. Both data show 
some signatures that are related to freezing and melt- 
ing of ^He in the nanopores. The freezing and melting 
occur at different temperatures and in finite tempera- 
ture ranges, unlike the first order transition of bulk ^He. 
We have identified P{T) at which ''He starts and termi- 
nates to freeze and thaw. The four data sets are plotted 
in Fig. [T] The reduction of freezing and melting tem- 
peratures observed above 3.7 MPa. Below 3.4 MPa, no 
signatures indicating freezing and melting were observed. 
The liquid-solid boundary below 0.8 K should therefore 
be located between 3.7 and 3.4 MPa, meaning that the 
freezing line is nearly flat and the entropy difference be- 
tween the solid and nonsuperfluid phases is small. This 
fact strongly suggests that the nonsuperfluid phase is a 
sort of an ordered state. 

To clarify the nature of the low-entropy nonsuper- 
fluid state, we have conducted the heat capacity 
measuremenl^'^. In Fig. [3ja) is shown the heat capacity 
data, in which ''He in the nanopores is liquid. A broad, 
but substantial peak is found in each heat capacity. The 
peak temperature Tb indicated by arrows and the peak 
height decrease as P increases. We plot Tb on the phase 
diagram of Fig. [TJ Obviously, the "Tb line" is located 
about 0.2K below the A line, and is parallel to the A line. 

The heat capacity peak is reminiscent of the super- 
fluid size effect in ''He in various restricted geometrieai^. 
However, the system exhibits no superfluid transition at 
and just below Tb. This is clearly seen in Fig. [Sjb), 
the data of the frequency shift in the torsional oscillator 
measurement. In Fig. [Htb), small upturns seen in both 
data around 2 K are due to the superfluid transition of 
the bulk liquid in the open space of the cell. The large, 
abrupt increase at lower temperatures indicates the su- 
perfluid transition of ^He confined in the nanopores. The 
superfluid transition temperature Tc is much lower than 
Tb, and it decreases progressively with increasing pres- 
sure. The remarkable difference in the behaviors of two 
characteristic temperatures is obviously seen in Fig. [1] 



B. Superfluid Properties 

Torsional oscillator technique'^'^ is based on a simple 
relationship that the frequency shift A/ is proportional 
to the superfiuid density p^. Therefore, A/(T) should 
contain essential information for understanding the na- 



ture of superfiuidity. In the next section we focus on the 
behavior of Ps(T). Here we mention some features in the 
A/(T) curves in two density regions: (1) adsorbed films 
to fiUed-pore states, (2) pressurized states at < P < 3.4 
MPa. 



1. Film States 

The adsorbed film shows the superfluid response 
when the coverage exceeds the critical coverage Uc = 
19.9 /xmol/m^. The superfluid transition temperature 
Tc increases almost linearly with the superfluid cover- 
age n — ric, and has a maximum at ninii — 33 /imol/m^, 
at which the pore is filled with ''He. It should be noted 
that the amount of the nonsuperfluid (i.e. "inert") layers 
adjacent to the pore walls are larger than that of super- 
fluid liquid under ambient pressure. The "effective" pore 
diameter for the superfluid part is therefore reduced to 
about 1.5 nm. 

From the slope of f{n) in the nonsuperfluid and super- 
fluid states, we obtained the ratio of undetected super- 
fluid mass to total superfluid mass, the so called x factor, 
to be 0.1. This value is much smaller than 0.33 in the 
case of Vycor, but larger than the obtained value by MT 
for the similar Gelsil glass, 0.06. These results indicate 
that the nanopores in Gelsil are more tortuous than the 
pores in Vycor. 

In Fig. [21 we show the typical superfluid frequency 
shift normahzed by the shift at K, A/(T)/A/(0), which 
is equal to the superfluid fraction ps (T) / p, together with 
the similar data of ^He film in 2.5-nm Gelsil by MT and in 
Vycor. The A/ data set at various coverages are shown in 
the inset of Fig. [IK a). There exists substantial difference 
in the temperature dependence of ps{T)/p among three 
experiments. Our Gelsil data lies between the Vycor data 
and Gelsil data by MT, and possess the features of these 
two systems. In the Vycor case, Ps{T) is proportional 
to T^ at low temperatures, and show a bulk-like critical 
behavior near Tj^. The T^ behavior suggests that one- 
dimensional phonons are the dominant low-energy exci- 
tations. In the Gelsil experiment by MT~, ps{T)/p is 
also fitted to T^ at low temperatures at T < 0.4 - 0.8Tc. 
As shown in Fig. [D(b), also in our Gelsil the normal fluid 
density obeys approximately T^ law at low temperature 
regions, as shown in the inset. Near Tc, our Ps{T) resem- 
bles the Vycor data, although the data are not enough 
to accurately determine the critical exponent. 



2. Liquid under Pressure 

Next we mention the results obtained in the pressur- 
ized states, where liquid ''He fills the nanopores^. In Fig. 
[H^b) and (c) are shown the A/(T) data. The data in 
Fig. [2Ib) are obtained by subtraction of the empty cell 
background at bulk superfluid trantision temperatures. 
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FIG. 2: Torsional oscillator results in the adsorbed film states, 
(a) Comparison of the normalized superfluid fraction {ps/p 
versus T/Tc) to the data of *He in other 2.5-nm Gelsil^" and 
in Vycor^^. The adopted data have similar Tc. Inset: Fre- 
quency shift A/(r) for eight coverages, (b) Log-log plots of 
the normal fluid fraction pn{T) / p for four coverages. The solid 
lines are the best powerlaw fittings. Inset: Powerlaw expo- 
nents a obtained in the fitting pn{T)/p = aT°' as a function 
of n. 
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We have found that also in pressurized liquid A/(T) 
obeys powerlaw at low temperatures. At P < 1.7 MPa 
the normal fluid fraction Pn/ p is best fitted by Pn/p oc 
with the exponent /3 ranging from 2.3 to 2.5. At higher 
pressures, a T-linear behavior emerges. In order to see 
the crossover from the nearly parabolic to linear temper- 
ature dependence, we have fitted the normal fraction to 
the sum of T-linear and square terms, pn/p — aT + bT^. 
The obtained coefficients a and b are plotted in Fig. |4l 
Obviously, the T-linear term dominates Pn/p above 2.3 
MPa. We will discuss the origin of the T-linear term in 
the next section. 



FIG. 3: (a) Heat capacity of liquid *IIe in the Gelsil nanopores 
for eight pressures. Arrows indicate the peak temperatures Tb 
that are interpreted as the LBEC formation temperatures. 
See recent publicationAiS for method of extraction of the heat 
capacity in the nanopores. (b) Torsional oscillator frequency 
shifts A/(r) at pressures below the bulk melting pressure. 
The shifts starting around 2 K are the contribution from 
bulk liquid in the sample cell, (c) Af{T) at pressures above 
the bulk melting pressure. The n-shaped anomalies are anti- 
crossing resonances resulting from the coupling to superfiuid 
fourth sound. 
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FIG. 4: The coefRcients a and b obtained from the hnear- 
parabola fitting A/(0) - A/(r) = aT + bT"^ for the data 
taken at nine pressures. 



III. DISCUSSION: LOCALIZED BEC AND 
PHASE FLUCTUATION 

A. The Localized BEC 




FIG. 5: A cartoon showing the formation of locaUzed BECs 
(LBEC) (shown as white) in a porous glass substrate (dark 
grey). ^He atoms form many small EEC's at the wider regions 
(especially intersections of the pores), where the atoms can ex- 
change actively. The phase of each BEC is illustrated by thick 
arrows. Since no phase coherence exist among the LBECs due 
to the hard core of *He suppressing the spatial exchanges at 
the narrower regions (light grey), the whole system exhibits 
no superfluidity on macroscopic length scale. As temperature 
is lowered, thermal phase fluctuations are diminished, and the 
system should undergo a macroscopic superfluid transition at 
some temperature Tc. 



We have proposed in the previous pubUcations that 
the QPT behavior, i.e. the anomalous reduction in su- 
perfluid Tc , results from the locahzation of Bose-Einstein 
condensates in the nonsuperfluid state^i^'S. We believe 
that this conjecture is now proven by the heat capacity 
measurement^!^. 

The idea of the localized BEC (LBEC) is shown in a 
cartoon of Fig. [5] When liquid ''He is confined in the 
nanopores, the BEC transition temperature should be 
reduced below bulk T\ due to the size effect. Around 
a certain temperature below many BECs grow from 
large pores or intersections of pores, in which ^He atoms 
can exchange frequently their positions. The size of the 
BECs is roughly limited to the pore size. The atom ex- 
change between the BECs via the narrow regions of the 
pores are interrupted, because ^He atom has a hard core. 
For the movement of one ^He atom, the surrounding ^He 
atoms act as a potential. The lack of the atom exchanges 
causes fluctuations in phase of the superfluid order pa- 
rameter. Therefore, no phase coherence exists among the 
BECs, and the whole system has also no global phase co- 
herence and does not exhibit superfluidity that can be 
detected by macroscopic and dynamical measurements 
such as torsional oscillator. As the temperature is fur- 
ther decreased, the phase coherence between the localized 
BECs grows, and macroscopic superfluidity, which is de- 
tected by torsional oscillator technique, is realized when 
most of the BECs coalesce. 

The heat capacity peak provides a definite evidence for 
the formation of LBECs at Tb • Broad peak structure in 
heat capacity is a common feature of ^He in restricted 
geometries, and was recognized as a manifestation of su- 



perfluidity and BEC. The temperature dependence of the 
heat capacity (the shape of the peak) of ''He in Gelsil 
agrees semi-quantitatively with that in restricted geome- 
tries such as Vycor. In our ^He-Gelsil system, however, 
the superfluid is much lower than the peak tempera- 
ture Tb j so the nanoscale BEC occurs around Tb without 
macroscopic superfluid transition. 

Heat capacity peak without macroscopic superfluid- 
ity has been observed in liquid "^Hc droplets formed in 
metal foils^^. In this case each droplet that is several 
nanometers in diameter is perfectly independent, and the 
droplets never exhibit superfluidity in macroscopic sense. 
The situation of ^He in nanoporous Gelsil is rather sim- 
ilar to this droplet system. 

In the abovementioned LBEC scenario, the smallness 
of the pore size is only essential to the QPT behavior. 
It has been pointed out that disorder or randomness 
in porous structures produces Boson localization called 
Bose glass state^'^. In the Bose glass state, the conden- 
sates localize at the local minima of the random poten- 
tials, and macroscopic phase coherence is lost by the lo- 
calization of atoms as in the case of narrowness-induced 
LBEC. Kobayashi and Tsubota have recently studied su- 
perfluidity of "^He confined in a 3D random model poten- 
tial taking account of the feature of our 2.5-nm GelsiU^. 
They found that superfluidity disappears above 4.2 MPa 
due to the localization of the BECs. It is in close agree- 
ment with our observation. 
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B. Effects of Phase Fluctuations 

In the LBEC state, phase of the superfluid order pa- 
rameter in each LBEC is fluctuating. This phase fluctu- 
ation can contribute to the superfluid properties below 
Tc. We propose that the T-linear behavior in the super- 
(or normal) fluid density observed at high pressures (Fig. 
2]) is the manifestation of the phase fluctuations that are 
induced thermally (classically). 

The effects of phase fluctuation has been studied in 
Josephson junction arraysi^ and granular metal film^i^, 
which show a superconductor - insulator quantum phase 
transition by controlling some experimental parameters 
such as magnetic field. It has also been proposed in 
the field of high-Tc cuprates that the phase fluctuations 
play an important role on the properties of underdoped 
regimes. Emery and Kivelson (EK)^^ argued that low 
carrier-density superconductors such as high-Tc cuprates 
are characterized by a small phase stiffness, and conse- 
quently the large phase fluctuations dominate notably 
the superconducting properties of underdoped regimes. 
The emergence of the "pseudo-gap" states is caused by 
the local Cooper pairing without global phase coher- 
ence throughout the sample, and the onset of long range 
phase order determines the true superconducting transi- 
tion that is detected by macroscopic means. This pro- 
posed mechanism is exactly the same as the LBEC pic- 
ture in the ^He-nanopore system. The LBEC state just 
corresponds to the EK pseudogap state. 

The superfluid systems that are controlled by phase 
fluctuations possess the following characteristics22i2^: 

1. The superfluid density ps is low, i.e. the phase stiff- 
ness (helicity modulus) is small even at K. 

2. The local order occurs at higher temperature than 
the long-range phase ordering. 

3. If the phase fluctuation is thermally excited, ps is 
proportional to T at low T . 

4. The long-range ordering is proportional to ps- 

In high-Tc cuprates, the T-linear behavior was ob- 
served by the measurement of penetration depth^SiSi. Al- 
though it is also attributed to the d-wave nature of the 
gap function, XY models with classical phase fluctuations 
reproduce quite well the overall temperature dependence 
of p;™. The smallness of and the proportionality be- 
tween Ps and Tc was also conflrmed and stressed as an 
important characteristic of various exotic superconduc- 
tors by Uemura et a L'^ '- i —: The proposed "universal" 
relation between Tc and muon relaxation rate was later 
reinterpreted by EK as an upper bound of Tc given by 
the phase-order temperaturei^. 

All the abovementioned features of the phase- 
fluctuation model are actually observed in the "^He-Gelsil 
system we studied. As is shown in Figs. [3] and [H the 
T-linear behavior in ps becomes prominent at pressures 
higher than 2.3 MPa. This behavior strongly suggests the 
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FIG. 6: Superfluid transition temperature Tc of pressurized 
liquid as a function of the frequency shift A/ at 10 mK. 
This plot corresponds to the Uemura plot for unconventional 
superconductora^ii^. 



existence of classical phase fluctuations which dominates 
the normal fluid component. It is also noted that the 
overall shape of the A/(T) curve bears striking resem- 
blance to the superfluid density of measured in cuprates 
and the calculated one in the 3D XY modeP^. 

Moreover, a plot of Tc versus ps (the so-called Ue- 
mura plot^iiS^) in Fig. [6] clearly have tendencies that 
at T > 2.3 MPa ps becomes small and approximately 
proportional to Tc. The emergence of T-linear term in ps 
correlates to the change in the slope of Tc — ps curve. 

The accuracy in the determination of T-linear coeffi- 
cient and Tc — ps relation near Pc are degraded in our 
current torsional oscillator measurement because of the 
small Ps (small signal-to-noise ratio) and the coupling of 
oscillation to fourth sound. Measurements of ps by other 
techniques such as fourth sound resonance method will 
be essential. 

The idea of LBEC gives a new perspective to a number 
of experimental studies of ^He in restricted geometries. 
The detailed torsional oscillator and specific heat stud- 
ies by Reppy and coworkers2SiS2i2^ shows that the su- 
perfluid transition occurs at slightly lower temperature 
than the temperature of the broad speciflc heat peak. 
At superfluid Tc, an extremely small peak is additionally 
observed. As well as in the Gelsil case, the broad peak 
is attributed to the formation of LBECs and the macro- 
scopic superfluid transition occur at Tc. The LBEC pic- 
ture in the *He-Vycor system is also supported by the 
neutron and ultrasound experiment conducted by Clyde, 
Mulders and coworkers, in which the roton signals are 
observed above Tc determined by ultrasound. Thus, the 
" separation" of BEC and superfluid transition should be 
a universal characteristic of ^He in nanoporous media. 

The T-linear superfluid density was observed in ^He 
filled in packed powders^ and in 2.5-nm CelsiliS at am- 
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bient pressures. In these studies the T-hnear behavior 
was attributed to the effect of "zero-dimensional" (OD) 
phonons. However, the OD phonons in the 3D connected 
nanopores are hard to imagine. It is rather reasonable 
to interpret as the effect of phase fluctuations. Then a 
question arises: Why is the T-linear behavior much more 
prominent in packed powder and MT's Gelsil than ours? 
It is conjectured that difference in pore structure influ- 
ence the Boson localization and thus the phase fluctua- 
tions. Further studies including detailed characterization 
of porous materials are obviously intriguing. 

IV. SUMMARY 

In summary, we have determined the anomalous phase 
diagram of "^He confined in the 2.5-nm Gelsil nanopores. 
It is ultimately proven by torsional oscillator and heat 



capacity studies that BEG and superfluidity take place 
at separate temperatures. Key physics to understand the 
phase diagram and superfluid properties is localization 
of Bose-Einstein condensates caused by confinement or 
disorder. Striking similarity to the superfluid behavior 
in high-Tc cuprates may also be an important clue to 
elucidate the mechanism of quantum phase transition. 
The entire phase diagram will be understood basically in 
terms of the phase fluctuation model that was proposed 
by Emery and Kivelson in the interpretation of phase 
diagrams of high-Tc cuprates. 

Our study shows that ^He in nanoporous media is an 
illustrative example of strongly correlated Bosons in po- 
tential, which produce intriguing quantum phenomena. 
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